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^S’his  report  documents  the  development  of  a  three-dimensional,  elastic- 
plastic,  work-hardening  constitutive  model  based  on  the  theory  of  effective 
(intergranular)  stress  for  the  purpose  of  simulating  the  response  of  isotropic 
three-phase  (partially  saturated)  earth  materials.  This  model  can  be  used  to 
perform  effective  stress  ground  shock  analyses. 
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20.  ABSTRACT  (Continued). 

— The  model  is  actually  a  marriage  ol  two  single-phaBe  elastic-plastic  models 
which  simulate,  respectively,  the  drained  and  undrained  stress-strain  behavior 
of  earth  media,  and  whose  coefficients  can  all  be  determined  experimentally 
using  conventional  laboratory  testing  techniques.  Fitting  the  model,  therefore, 
requires  an  extra  set  of  laboratory  *.  'nts  (the  drained  tests)  over  that  required 
to  fit  a  contemporary  (total  stress)  ground  shock  model.  , 

For  undrained  applications  (where  only  the  applied  total  stresses  are 
usually  known),  the  model  predicts  total  stresses  and  ground  deformations  as 
well  as  the  effective  stresses  and  pore  pressures.  Pore  pressure  is  an  output 
from  the  model,  not  an  input.  Laboratory  pore  pressure  data  obtained  during 
conduct  of  the  undrained  testa,  however,  are  extremely  useful  for  verification 
of  the  model. 

The  model  is  examined  under  simulated  drained  and  undrained  laboratory  test 
conditions  to  depict  its  various  features  and  illustrate  the  techniques  for  fit¬ 
ting  its  functional  coefficients.  In  addition,  given  a  set  of  stress  increments 
or  a  set  of  strain  increments,  two  general  methods  are  presented  for  exercising 
the  model  to  compute  the  resultant  total  stresses,  effective  stresses,  pore 
pressures,  and  deformations. 

- — The  feasibility  of  using  the  new  model  for  effective  atreBB  ground  shock 
analyses  is  demonstrated  numerically  by  solving  a  simplified  tvo-dimenBional 
airslap  problem  for  a  saturated  sand  halfspace  using  an  effective  stress  version 
of  the  LAYER  code.  The  real  power  of  this  computational  tool  has  yet  to  be 
investigated  by  applications  against  actual  field  experiments. 

'■  >  Appendix  A  reviews  the  fundamentals  of  elastic-plastic  models. 
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CHAPTER  1 


INTRODUCTION 


1.1  BACKGROUND 

Earth  materials  are  multiphase  systems  (Reference  l)  that  consist 
in  general  of  solid  particles  (possibly  cemented) ,  water,  and  gas.  The 
intrinsic  response  of  such  materialr  to  externally  applied  loads  Is 
extremely  complicated.  To  model  this  response  for  a  particular  material, 
one  must  resort  to  the  theory  of  continuum  mechanics  and  have  available 
an  appropriate  constitutive  relation.  The  solution  of  earth  structure 
problems  then  becomes  a  mathematical  formalism  that  can  be  achieved 
numerically  or  by  other  means. 

In  recent  years,  primarily  under  the  sponsorship  of  the  Defense 
Nuclear  Agency  (DHA),  considerable  progress  has  been  made  in  the  devel¬ 
opment  of  mathematical  models  for  the  study  of  ground  shock  effects  in 
soil  and/or  rock  media  (References  2  through  12).  The  models  are  used 
in  two-dimensional  (2D)  computer  code  calculations  that  help  define  the 
ground  shock  environment  of  current  and  planned  hardened  defense  facili¬ 
ties.  The  requirement  for  more  advanced  models  can  be  attributed  to  the 
desire  for  a  better  mathematical  approximation  of  the  observed  stress- 
strain  properties  of  earth  materials.  Consequently,  several  quite 
complicated  nonlinear  elastic-ideally  plastic  (References  2  through  k), 
variable  modull-type  (References  5  through  8),  and  nonlinear  elastic- 
plastic  work-hardening  (References  9  through  12)  constitutive  models 
have  been  developed  and  used  by  the  DNA  soil  and  rock  mechanics  communi¬ 
ties.  A  detailed  summary  of  the  advantages  and  disadvantages  of  each 
type  of  model  is  given  in  Reference  3.  A  brief  historical  development 
of  these  models  is  given  in  Reference  12.  All  of  the  above  models  are 
three-dimensional,  either  isotropic  or  transverse-isotropic,  and  can 
simulate  to  various  degrees  the  observed  highly  nonlinear  and  hysteretic 
behavior  of  earth  media,  with  some  of  them  predicting  shear- induced 
volume  change;  however,  they  all  have  in  cowman  that  they  can  simulate 
only  single-phase  systems. 


Treatment  of  two-phase  systems  (water  and  solid)  for  fluid- 
saturated  earth  materials  within  the  framework  of  elastic-plastic 
constitutive  models  was  first  attempted  by  Roscoe  and  his  coworkers  at 
Cambridge  University  (Reference  13).  They  successfully  developed 
several  two-dimensional  models  for  describing  the  stress-strain-pore 
pressure  response  of  saturated  clay.  The  solid  skeleton  and  the  fluid 
constituent  of  their  two-phase  system  (saturated  clay)  were  assumed  to 
be  incompressible,  thus  allowing  the  treatment  of  undrained  loading 
conditions  by  specifying  ds^  =  0  (de^  =  increment  of  volumetric 
strain).^-  The  assumption  of  the  incompressibility  of  the  solid  skeleton 
and  the  fluid  constituent  is  reasonable  since  in  most  conventional 
geotechnical  problems  the  loading  conditions  are  static  and  the  pressure 
levels  of  interest  are  very  small.  In  1976  and  1977,  several  isotropic 
two-phase  constitutive  models  for  saturated  cohesionless  soils  were 
developed  (References  lU  and  15)  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES).  These  models  also  used  the  concept  of 
de^  *  0  to  simulate  undrained  test  conditions  for  fluid-saturated 
granular  materials.  In  contrast  to  the  Cambridge  models  (Reference  13), 
however,  the  WES  models  are  three-dimensional  and  provide  greater  flexi¬ 
bility  in  fitting  test  data  for  a  broad  range  of  behavior.  The  models 
documented  in  References  13,  lU,  and  15  provide  the  means  to  perform 
effective  stress  analyses  for  problems  involving  static-type  loading 
conditions.  Because  of  their  assumption  that  the  solid  skeleton  and  the 
water  are  incompressible  (i.e.,  dc^  *  0  ),  however,  they  are  not 
suitable  for  problems  involving  transient  or  dynamic-type  loading  con¬ 
ditions,  such  as  ground  shock  problems.  The  elimination  of  the  assump¬ 
tion  of  incompressibility  of  the  solid  skeleton  and  water  not  only  pro¬ 
vides  a  means  of  solving  dynamic  problems  but  also  provides  a  means  for 
simulating  partially  saturated  materials  (i.e.,  solid  skeleton,  water, 
and  gas).  The  model  developed  in  this  report  closely  parallels  the  devel¬ 
opments  reported  in  Reference  lU  with  the  exception  that  the  assumption 
of  incompressibility  of  the  solid  skeleton  and  water  is  eliminated. 

^  For  convenience,  symbols  and  abbreviations  are  listed  and  defined,  in 
the  Notation  (Appendix  B). 
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The  development  of  the  present  model  ia  based  on  the  fact  that  cur¬ 
rent  single-phase  and  tvo-phase  constitutive  models  vill  each  predict 
a  different  deformation  path  for  a  given  stress  loading ,  with  the  tvo- 
phase  model  achieving  the  better  agreement  with  physics  and  data.  By 
implication ,  a  two-  or  three-phase  constitutive  model  vill  also  predict 
better  ground  motion  histories  than  the  current  crop  of  single-phase 
ground  shock  models.  To  vhat  extent  and  in  vhat  situations  remain  to 
be  seen. 

1.2  OBJECTIVE 

The  overall  objective  of  this  study  vqb  to  develop  a  completely 
general ,  three-dimensional,  elastic-plastic,  vork-hardening  constitutive 
relationship  for  simulating  the  behavior  of  isotropic  three-phase  earth 
materials. 

1.3  SCOPE 

One  of  the  constraints  of  this  study  vsb  that  the  model  should  be 
in  a  form  suitable  for  use  vith  current  finite-difference  techniques 
(such  as  the  LAYER  code)  for  the  computation  of  ground  motions,  total 
and  effective  stresses,  and  pore  pressures  produced  vithin  earth  masses 
by  explosion-induced  ground  shock.  The  vork  involved  five  steps: 

(l)  extending  the  mathematical  model  documented  in  Reference  1**  to  allov 
for  volume  changes,  (2)  fitting  the  model  to  available  drained  and  un¬ 
drained  mechanical  property  test  data,  (3)  validating  the  model  by 
calculating  the  pore  pressure  responses  for  various  laboratory  test 
boundary  conditions  and  comparing  the  predictions  vith  actual  test  data, 
(M  devising  an  efficient  numerical  logic  for  calculating  time  histories 
of  effective  stress,  pore  pressure,  and  associated  ground  motions  for 
blast-oriented  problems,  and  (5)  incorporating  the  model  and  its  calcula- 
tional  logic  into  LAYER  and  performing  a  demonstration  calculation. 

The  concept  of  effective  stress  and  the  mechanical  behavior  of  soil 
are  presented  in  Chapter  2.  In  Chapter  3»  the  development  of  a  single¬ 
phase,  elastic-plastic  constitutive  model,  parallel  to  that  reported  in 
Reference  lU,  is  presented.  The  application  of  this  model  to  treat 
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multiphase  systems  1b  explained  in  Chapter  U .  The  quantitative  behavior 
of  the  multiphase  model  under  simulated  triaxial  test  conditions  is 
examined  in  Chapter  5.  Chapt-'r  b  describes  mathematical  fits  developed 
with  available  drained  and  undrainod  mechanical  property  test  data  for 
Campbell  Swamp  Band  and  presents  comparisons  of  model-predicted 
responses  with  the  measured  laboratory  behavior.  The  results  of  the  CP 
demonstration  calculation  are  presented  in  Chapter  ?.  Chapter  8 
summarises  key  aspects  of  the  model  and  offers  recommendations  for  its 
quantitative  evaluation.  Appendix  A  reviews  the  fundamental  basis  of 
elastic-plastic  constitutive  models  and  is  included  for  reference 
purposes  and  future  use. 
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CHAPTER  2 


EFFECTIVE  STRESS  CONCEPT  AND  MECHANICAL  BEHAVIOR  OF  SOIL 

Earth  material,  in  its  general  form,  is  composed  of  a  complex 
assemblage  of  discrete  particles  of  varying  shapes  and  orientation  in  a 
compact,  possibly  cemented,  array.  Theae  particles  may  range  in  magni¬ 
tude  from  the  microscopic  elements  of  a  clay  soil  to  the  macroscopic 
boulders  of  a  rock  fill.  The  voids  in  the  array  may  be  filled  with 
water  or  air  and  usually  contain  both.  Before  a  constitutive  model 
describing  the  behavior  of  these  materials  under  an  applied  stress 
system  can  be  developed,  it  is  necessary  to  consider  how  these  stresses 
are  distributed  among  the  several  components  comprising  the  aggregate 
and  to  understand,  in  general  terms,  the  mechanical  behavior  of  these 
assemblages.  The  emphasis  throughout  the  remainder  of  this  report  will 
be  on  earth  materials  that  are  better  described  as  soils  than  rocks; 
however,  the  model,  in  principle,  is  applicable  to  both. 

2.1  EFFECTIVE  STRESS  CONCEPT 

The  normal  stress  components  at  a  point  in  a  soil  body  may  be 
divided  into  two  parts  (Reference  16):  The  stress  carried  by  the  solid 
skeleton,  referred  to  as  the  effective  stress,  and  the  stress  carried 
by  the  pore  fluid,  referred  to  as  the  pore  pressure.  The  pore  pressure, 
in  turn,  can  be  divided  into  two  additional  parts:  the  stress  carried 
by  the  water  and  the  stress  carried  by  the  air.  According  to  Refer¬ 
ences  17,  lfl,  and  19,  total  stress^  can  be  expressed  (in  indicial 
notation)  as 

"u  '  *ij  *  [r. '  *(r»  -  r.)]5u  <?-u 

1  Indices  take  on  values  of  I,  2,  or  3.  A  repeated  index  is  to  be 
summed  over  its  range,  A  comma  between  subscripts  reprer-nts  a 
derivative.  Quantities  are  referred  to  rectangular  Carte*. an 
coordinates  . 
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where 


°1J  "  stress  tensor 
°iJ  ■  effective  stress  tensor 
Pa  *  Pore  air  pressure 

x '  oc^X'"  ^ 

Pw  *  P°r®  water  pressure 


to  the  pore  volume 


4ij  *  decker  delta  ■ 


!*•  *  - 

( 0 ,  i  i 


Equation  2.1  cw  be  rewritten 


°U  *  °ij  * 


(2.2) 


°r  th'  p°r'  *ir  >«■.»»  :nhr;r.:::r" th'  'r™ 

^  water  pressure,  i.e,, 


U  “  Pa  *  *  Pv) 


(2.3) 


Tor  a  fully  saturated  soil,  x  .  1  . 

X  -  0  .  ♦  *nd  for  a  completely  dry  soil, 

FOr  *  trU*U1  t«t  performed  on  .  , 

coordinate  ,„te.,.  th.  ^  ^ 


°.  "o’  ♦  u 

•  « 


°6  °r  "  °e  *  u  ■  ♦  u 


(2.M 


where  o  ,  0  . 

‘«*entlL  totL  .tr...°^nuTd°tlo"1Jr'  ‘I”  ^  *"d 

— effeetle,  co._tl*  ££  of 


lk 


saturated  and  partially  saturated  soils  tested  under  trlaxial  conditions 
is  discussed  in  the  following  section. 

2.2  MECHANICAL  BEHAVIOR  OF  SOIL 

The  mechanical  behavior  of  soils  subjected  to  externally  applied 
load,  is  quite  complicated.  Unlike  the  properties  of  most  engineering 
materials,  soil  stress-train  properties  are  greatly  affected  by  such 
factors  as  soil  structure,  degree  of  voids  saturation,  drainage  condi¬ 
tions  during  loading,  loading  rate,  loading  history,  and  current  stress 
state.  In  this  report,  the  tens  soil  structure  is  loosely  used  to 
recognise  that  "microvwlable."  such  ..  grain  size;  grain  site  distribu¬ 
tion;  grain  shape,  surface  texture,  and  mineralogy;  and  grain  orientation, 
packing,  and  cementation,  or  bonding,  all  play  an  interrelated  role  in 
the  complex  mechanical  response  of  soils.  For  practical  purposes  these 
micro variables  we  generally  not  Individually  characterized;  rather 
engineering"  ""lables  such  as  void  ratio  and  relative  density,  which 

*"  “"en“ble  to  “easurement  and  certainly  more  useful  to  analysts, 
and  which  account  in  a  lumped  sense  for  the  interplay  of  many  of  the 

micro  variables ,  are  used  to  help  characterise,  understand,  snd  predict 
the  soils  response. 

Figure  2.1  show,  a  typlcol  behavior  of  soil  subjected  to  a  hydro¬ 
static  state  of  stress.  It  is  clew  from  this  figure  that  soils,  i„ 
general,  exhibit  a  nonlinear  compacting  hydrostat. 

Figure  2.2  shows  a  typical  vwlety  of  stre.s-strain-pore  pressure 
response  curve,  manifested  by  saturated  soil,  tested  in  undrained  shew 
n  a  triaxial  compression  device.  The  three  specimens  were  first 
isotropically  consolidated  to  the  same  effective  mew  normal  stress 
evel  (point  2),  then  shewed  undrained.  He  shew  curves  mwked  ”2  - 
show  the  typical  response  of  a  normally  consolidated  clay  or  a  very 
loos,  sand.  The  curve,  marked  "2  -  5"  show  behavior  typical  of  an  over- 
consolidated  clay  or  a  very  dense  sand.  Within  the  extreme  limit,  of 


deformat  ion.  *U,t  -~.ur-.nt.  of  radial 


ijr 


these  loose  and  dense  soil  responses,  there  is  a  graduated  response, 
typified  herein  by  the  curves  marked  "2  -*•  U."  The  latter  response 
depends  on  the  state  of  compaction  (consolidation)  of  the  material.  It 
is  clear  from  this  figure  that  the  effective  stress  is  the  only  part  of 
the  total  stress  that  affects  soil  shear  strength. 

Figure  2.3  shows  typical  (qualitative)  stress-strain  response 
curves  for  soils  sheared  under  drained  triaxial  compression  conditions; 
i.e.,  the  curves  marked  "1”  represent  dense  sand  or  overconsolidated 
clay,  while  the  curves  marked  "2"  depict  response  typical  of  loose  sand 
or  normally  consolidated  clay. 

In  Chapter  3  the  mathematical  development  of  a  total-stress  (single¬ 
phase),  elastic-plastic,  isotropic  constitutive  relationship  that 
can  qualitatively  describe  most  of  the  behavior  shown  in  Figures  2.1 
through  2.3  is  presented.  This  model  is  unable  to  simulate  the  postpeak 
work-softening  shear  response  exhibited,  for  example,  by  material  1  in 
Figure  2.3. 


1  Typical  behavior  of  a  dry  or  drained  soil 
under  hydrostatic  loading  and  unloading. 


IT 
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Figure  2.3  Typical  behavior  of  saturated  soils 
tested  under  drained  tri axial  test 
conditions . 
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CHAPTER  3 

ELASTIC-PLASTIC  CONSTITUTIVE  MODEL 

The  model  described  in  this  chapter  is  a  fundamental  element  of  the 
multiphase  model  discussed  in  Chapter  fc.  Hence,  the  understanding  of 
this  model  will  smooth  the  transition  to  the  multiphase  model. 

The  basic  theory  of  elastic-plastic  constitutive  models  for  single- 
phase  (solid)  materials  is  presented  in  detail  in  Appendix  A.  The 
elastic  behavior  of  these  models  is  defined  by  Equations  A. 12  through 
A.16  of  Appendix  A.  The  plastic  behavior  is  described  by  Equations  A. 17 
through  A. 28.  The  complete  elastic-plastic  description  is  expressed  by 
Equation  A. 29  and/or  Equation  A. 30.  These  equations  are  used  herein, 
with  selected  mathematical  forms  of  the  various  response  functions  con¬ 
tained  in  the  model,  to  describe  how  most  of  the  typical  soil  responses 
presented  in  Chapter  2  can  be  simulated  with  a  more  or  less  conventional 
elastic-plastic  model. 

3.1  ELASTIC  BEHAVIOR 

The  behavior  of  the  model  in  the  elastic  (recoverable)  range  is 
governed  by  the  elastic  bulk  and  shear  moduli  (Equation  A. '6  of  Appen¬ 
dix  A),  The  elastic  bulk  modulus  describes  the  unloading  stress-strain 
response  of  a  hydrostatic  compression  test  (Figure  3.1).  It  is  sug¬ 
gested  that  for  most  isotropic  earth  materials  the  elastic  bulk  modulus 
can  be  taken  as  a  function  of  the  mean  normal  stress,*  P  ,  or  the  first 
invariant  of  the  stress  tensor"  (J^  ■  3P).  The  following  expression 
is  chosen: 


K 


1  -  K, 


Ki  "'  (-Vi)] 


(3.1) 


The  elastic  bulk  modulus  could  also  be  a  function  of  the  plastic 
volumetric  strain. 

In  this  report  tension  is  considered  as  negative. 
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where 

*  initial  elastic  bulk  modulus 
and  *  material  constants 

Equation  3.1  is  illustrated  in  Figure  3.2.  Equation  3.1  together  with 
Equation  A. 13  of  Appendix  A  indicates  that  the  material  constants 
Kj  ,  ,  and  K0  can  be  readily  determined  experimentally  from  the 

unloading  hydrostatic  compression  test  results*  as  illustrated  in 
Figure  3.1. 

3 

The  elastic  shear  modulus  must  account  for  the  curvature  observed 

ii.  unloading  stress  difference-strain  difference  results  obtained  from 

triaxial  compression  tests  (Figure  3.3).  For  this  report,  the  elastic 

shear  modulus  is  assumed  to  be  a  function  of  the  second  invariant  of  the 

stress  deviation  tensor,  J0  ,  and  the  plastic  volumetric  strain, 

P 


TTcT  [1  ’  G1  «xp(-g?^2>3  ♦  °3  E1  -  ekk)]  (3-2) 


where 


■  initial  elastic  shear  modulus 


G1  *  Gi  *  G3  *  anc*  Gj1  "  material  constants 


Equation  3.2  is  illustrated  in  Figure  3.1*.  Equation  3.2  together  with 
Equation  A. 13  of  Appendix  A  indicates  that  the  material  constants 
Gi  ’  Gi  ’  G->  »  Gj  *  ant*  ^  can  P*  readily  determined  from  the 

slopes  of  experimental  unloading  stress-strain  curves  obtained  from  a 
series  of  triaxial  shear  tests  conducted  at  different  confining  pres¬ 
sures  (Figure  3.3). 


3  The  functional  forms  of  the  bulk  and  shear  moduli  (Equations  3.1  and 
3.2)  could  include  ax?re  terms  and,  hence,  provide  more  flexibility 
in  fitting  the  behavior  of  any  specific  material  (ef.  Reference  10). 
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3.2  PLASTIC  BEHAVIOR 

or  Jir.rr  b‘havior' ,he  io*dine  i  (— »!« *.  ^ 

l  “ to  be  iaotropic  “d  *°  °°“iat  °f  «» 

-5).  an  ultimate  failure  envelope  which  eervee  to  limit  the 

T7  “tlln,b1'  b*  th«  and  an  elliptical!* 

P«h  strain-hardening  yield  surface  that  produces  plastic  volumetric 

ZtTlr  Tl’Z  U  "°VM  <R'ferenCe  10K  "•  envelope 

£  ‘  D  fU”0tIOn  ls  “SUJ'd  t0  *  °f  «•  Seneralited 

Prager-Drucker  type  and  is  mathematically  described  by 


f,Jr  ■G?  *  ^  -  U  -  c  e*p  (-BJ^)J  .  o  (3 

•ni  the  strain-hardening  yield  surface  is  described  by 

F<V  <)  -  [Jx  -  LU)]2  .  -  [XU)  .  lU)]2  .  0  (j 

“fcore  A  .  B  ,  and  C  are  material  constants  (Figure  3.5);  R 

intersectiWhiChfWi11  “  d'fI"'d  ‘*l0Vi  ^  LU)  *«»  «* 

rsectlon.  of  the  hardening  surf.ee  with  the  J.  „is  and  the 

failure  envelope  ffj  /r  \ 

parameter  wht7  *’  2  '  “1  <  is  the  hardening 

parameter,  which  generally  is  a  function  of  the  history  of  plastic 
volumetric  strain,  £  .  For  _t  soil>>  .  _  ^  ^  ^ 


(3.3) 


(3.M 


«  ■  Ca 


(3.5) 


Eguation  3.5  allow,  the  elliptic  hardening  surface  to  eap»d  and 

ana. ”'»>!*  *'  *°  tr“*Ut"  r'1,tlYr  t0  'he  origin  of  the  J  ,  /J 

Z  *  '  th'  harde”lae  ”arf*"  <«**»  3.5)  V..  chosen  aoW 
«*•«  at  it.  intersection  with  the  failure  envelope  i.  horisontal. 


■ 


This  condition  is  guaranteed  by  the  following  relationships  between  <  , 
L(«c)  ,  and  x(k)  :* 


!t(ic)  if  1(k)  >  0 

(3.6) 

0  if  1(k )  <_  0 

*  *  ekk  *  W{1  "  exp  ♦  W1[x(k) 32  exp  [-D^U)]  (3.7) 

X(<)  *  t(<)  +  R  { A  -  C  exp  [-Bt(<)]}  (3.8) 

where  D  ,  1)^  ,  and  are  material  constants,  and  W  ic  also  a 

material  constant  which  defines  the  maximum  plastic  volumetric  compac¬ 
tion  that  the  material  can  experience  under  hydrostatic  loading  (Figure 
3.1).  The  parameter  R  ,  in  Equations  3.1*  and  3.8,  is  the  ratio  of  the 
major  to  the  minor  axes  of  the  elliptic  yield  surface  (Figure  3*5).  The 
value  of  R  depends  on  the  state  of  compaction  of  the  material.  For  a 
contractive  material  (i.e.,  loose  sand  or  normally  consolidated  clay, 
curves  marked  "2  *  3"  in  Figure  2.2),  the  value  of  R  is  greater  than 
1/a  where  a  ■  CB  exF  (-BJ^  is  the  slope  of  the  failure  envelope 
(Equation  3.3).  For  a  dilative  material  (i.e.,  dense  sand  or  overcon¬ 
solidated  clay,  curves  marked  •*  5"  in  Figure  2.2),  the  value  of  R 
is  less  than  l/o  .  R  *  i/a  corresponds  to  the  curves  marked  "2  ■*  1»M 
in  Figure  2.2.  These  variations  in  the  parameter  R  can  be  accounted 
for  by  the  following  functional  relation: 

R.  2 

R  ■  p  +  p  -  (  1  [-R0L(k) 3)  ♦  R^  exp  {-RjjLU)  -  R,.]  )  (3.9) 


*  The  mathematical  form  of  Equation  3.7  depends  on  the  specific  mate¬ 
rial  being  modeled.  The  author  believes,  however,  that  the  form 
presented  by  Equation  3*7  is  suitable  for  modeling  most  soils. 
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where  ,  R1  ,  »  R^  •  R^  ,  and  R^  are  material  constants  that 

can  be  determined  by  a  trial-and-error  process  of  fitting  the  model  to  a 
variety  of  laboratory  stress-strain  data. 

3.3  SUMMARY 

In  summary,  there  are  five  potential  functions  (two  elastic  and 
three  plastic)  that  describe  the  complete  behavior  of  the  single-phase 
model.  These  functions  are  summarized  in  Table  3.1.  In  addition,  21 
material  constants  are  used  to  characterize  these  functions  in  the 
present  model.  They  too  ave  summarized  in  Table  3.1. 

In  the  next  chapter,  it  is  shown  how  the  present  single-phase  model 
can  be  adapted  to  simulate  multiphase  soil  materials. 
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Figure  3.5  Proposed  yield  surfaces  for  the  elaatic-plnstic 
strain-hardening  model. 


CHAPTER  1» 


THE  TREATMENT  OF  A  MULTIPHASE  SYSTEM 

The  three-phase  soil  model  is  patterned  on  the  concepts  of  effec¬ 
tive  stress  presented  in  Section  2.1  and  is  really  quite  simple. 

First  of  all,  two  single-phase  models  of  the  type  developed  in  Chapter  3 
are  required.  One  of  these  is  used  to  characterize  the  effective  stress 
response  of  the  soil  of  interest,  and  the  other  is  used  to  characterize 
the  toted  stress  response,  as  determined  by  appropriate  laboratory 
tests.  The  only  remaining  step  is  to  combine  these  two  models.  This 
is  accomplished  by  imposing  two  constraints:  (l)  compatibility  of 
toted  volumetric  strain  increments  emd  (2)  compatibility  of  individual 
effective  and  total  deviatoric  stress  increments  (dSj^  -  dS^).  The 
first  constraint  leads  to 

J  ■  Jj  +  u  (U.l) 


where  is  the  first  invariant  of  the  effective  stress  tensor;  and 

the  second  leads  to 


•ij  *  "'u 


which  1 8  Equation  2.2.  The  resulting  model  is  thus  able  to  predict 
deformations,  total  and  effective  stresses,  and  pore  pressures  for  real 
three-phase  media.  The  predicted  pore  pressures,  however,  are  the 
overall  pore  pressures,  u  ;  the  model  dees  not  distinguish  between  the 
pore  air  pressure  P  and  the  pore  water  pressure  P  discussed  in 

^  R  V 

Chapter  2.  Consequently,  it  is  only  a  three-phase  model  in  a  "phenome¬ 
nological"  sense.  Perhaps  it  should  be  thought  of  as  a  "pseudo"  three- 
phase  model.  Be  that  as  it  may,  the  model  does  predict  observed 

1  The  pore  air  preesure  PH  and  the  pore  water  pressure  Pw  could  be 
calculated  from  Equation  2.1  if  an  appropriate  laboratory  test  is 
conducted  to  determine  the  value  of  x  • 


three-phase  material  response;  and  it  does  this  quite  well,  as  will  be 
demonstrated  in  the  next  chapter. 

Technically  speaking  then,  modeling  the  behavior  of  multiphase 
soil  systems  requires  two  separate  sets  of  material  constants,  such  as 
those  shown  in  Table  3.1.  The  first  set  must  reflect  the  effective 
stress  properties  of  the  soil  (i.e.,  the  properties  of  the  soil  skeleton 
alone)  and  must  be  determined  by  fitting  test  data  obtained  for  the 
material  under  fully  drained  conditions.  The  second  set  must  reflect 
the  total  stress  properties  of  the  soil  (i.e.,  those  of  the  skeleton- 
water-air  mixture),  which  must  be  determined  by  fitting  test  data 
obtained  for  the  material  under  completely  undrained  conditions.  The 
resulting  two  sets  of  mod«l  parameters  are  summarized  in  Table  .  1 . ^ 

Combining  these  two  models  through  volumetric  strain  compatibility 
and  effective  and  total  deviatoric  stress  compatibility  allows  for 
calculation  of  the  pore  pressure  and  deformation  response  and  the  total 
and  effective  stress  response  of  a  multiphase  system  subjected  to  given 
stress  or  strain  increments.  Either  one  of  the  following  two  proceduree 
can  be  used  for  this  purpose: 

1.  If  stress  increments  are  given, 

a.  Calculate  the  undrained  volumetric  strain  using  the 
second  set  of  response  functions  and  material  con¬ 
stants  listed  in  Table  U.l,  i.e.,  the  undrained 
model  parameters. 

h.  Impose  this  volumetric  strain  on  the  drained 
model  (i.e.,  the  first  set  of  response  functions 
and  material  constants  listed  in  Table  b.l)  and 
calculate  the  resulting  stress  path  and  associated 
material  response.  This  stress  path  is  the  effec¬ 
tive  stress  path  that  the  material  will  experience 
during  this  undrained  load  application.  The  pore 
pressure  is  simply  the  difference  between  the  total 
and  the  effective  normal  stresses. 

This  procedure  is  illustrated  in  the  following  diagram: 

2  Rote  that  1*2  material  constants  are  required  to  fully  define  the 
general  model  presented  in  Table  l».l.  Depending  on  the  actual  material, 
however,  fewer  than  1*2  constants  may  be  needed. 
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2.  If  strain  increments  are  given  (the  usual  case  in  code  calcula¬ 
tions), 

a.  Calculate  the  effective  stress  using  the  first  (drained) 
set  of  response  functions  and  material  constants  listed 
in  Table  U.l. 

b.  Calculate  the  total  stresses  using  the  second  set  of 
response  functions  and  material  constants  listed  in 
Table  U.l,  The  pore  pressure  during  this  undrained 
load  application  is  simply  the  difference  between  the 
total  and  the  effective  normal  stresses. 

The  following  diagram  illustrates  this  procedure: 


The  response  of  this  multiphase  material  model  subjected  to 
undrained  standard  triaxial  test  conditions  is  examined  in  detail  in  the 
next  chapter  using  the  first  of  these  procedures.  This 


exposition  of  the  model  will  help  the  reader  to  appreciate  both  its 
relative  simplicity  and  its  power. 


CHAPTER  5 


BEHAVIOR  OF  THE  MULTIPHASE  CONSTITUTIVE  MODEL  UNDER 
TRIAXIAL  TEST  CONDITIONS 


5.1  INTRODUCTION 

The  ability  of  the  new  model  to  simulate  the  response  of  multiphase 
soil  systems  can  be  more  clearly  understood  if  the  model  is  examined 
under  familiar  laboratory  test  boundary  conditions.  Since  most  of  the 
mechanical  testing  of  soils  for  engineering  purposes  is  performed  with 
the  triaxial  test  (TX)  apparatus,  it  is  appropriate  to  investigate  the 
model  under  both  drained  and  undrained  TX  conditions.  Adopting  the  z- 
axis  of  a  cylindrical  coordinate  system  (z,  r,  and  8)  as  the  axis  of  sym¬ 
metry  of  the  test  sample,  the  total  and  effective  stress  tensors  and  the 
total  strain  tensor  associated  vith  this  configuration  become, 
respectively 


where 

tively. 


(5.1) 


(5.2) 


(5.3) 


c  and  c  are  the  total  vertical  and  radial  strains,  re spec- 
i  r 

Equations  5.1  and  5*2  imply  homogeneity  of  stress  (oA  ■  o  ), 

o  r 


3T 


and  Equation  5.3  implies  homogeneity  of  deformation  (total  tangential 

strain  e  •  e  ).  The  variables  P*  ■  J* /3  (effective  mean  normal 
“  r  x  __ 

stress),  P  «  J^/3  (total  mean  normal  stress),  (the  second  invari¬ 

ant  oi'  the  effective  stress  deviation  tensor),  J.  (the  second  invari- 

d 

ant  of  the  total  stress  deviation  tensor),  and  e^/3  (mean  volumetric 
strain)  associated  vith  the  above  stress  and  strain  tensors  take  the 
following  forms 


o’  +  2o' 


p*  .  j./3  .  -L- - L 

(5.M 

a  ♦  2a 

P  -  Jx/3  -  - - 

(5.5) 

(o'  -  0»)2  (o  -  cr  )2 

j»  «  j*  »  — 3- — -J— >  —3 - 1 — 

2  2  3  3 

(5.6) 

1  AV  ekk  cz  +  2er 

3  V„  3  3 
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where  AV/V  is  the  volumetric  strain.  The  TX  test  is  generally  con- 
o 

ducted  in  two  phases:  the  hydrostatic  phase  and  the  shear  phase.  Both 
phases  can  be  conducted  either  drained  or  undrained.  They  are  discussed 
below. 

5.2  HYDROSTATIC  PHASE 
5.2.1  Drained  Condition 

During  the  drained  hydrostatic  phase  of  a  TX  test,  the  pore  pres¬ 
sure  is  always  zero,  and  the  behavior  of  the  soil  skeleton  alone  (i.e., 
the  effective  stress  behavior)  is  examined.  The  following  conditions 
are  obtained: 
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The  relation  between  the  increment  of  the  first  invariant  of  effective 
stress  and  the  elastic  volumetric  strain  increment  is  given  as  u«e 
Equation  A. 13  and  Set  1  of  Table  U . 1 ) 
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where  the  response  function  K  is  given  by  Equation  3.1,  Equation  3.1 

s 

is  substituted  into  Equation  5*10*  and  the  resulting  expression  la 
integrated  to  provide  the  following  relation  between  the  elastic 

JJ 

volumetric  strain  and  ,T* 
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The  relation  between  the  plastic  volumetric  strain,  i  ^  and  d '  is 
given  by  Equation  3.T*  where  *  for  thin  phase  of  the  test,  Is  r|^ 
and  x ( k )  is  ,  thus! 
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In  view  of  Equations  5*11  and  5*12,  the  total  volumetric  stratn  taken 
the  following  form 
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Equations  5.9  and  5.11  through  5.13  provide  a  complete  specification  for 
the  deformation  response  of  the  soil  skeleton  subjected  to  a  drained 
hydrostatic  test  (i.e,,  isotropic  consolidation). 


The  qualitative  behavior  of  the  model  during  a  drained  hydrostatic 
test  is  shown  in  Figure  5.1.  The  slope  of  the  curve  during 

virgin  loading  can  be  obtained  from  Equation  5*13  with  the  help  of 
Equations  5*12  and  3.1: 
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and  K  and  K  are,  respectively,  the  elastic  (Equation  3.1  and 
s  s 

Table  l*.l)  and  the  apparent  bulk  moduli  of  the  soil  skeleton  under 
drained  hydrostatic  loading.  The  second  and  third  terras  in  the  denomi¬ 
nator  of  Equation  5-1^  produce  a  softening  of  the  apparent  bulk  modulus 
due  to  plastic  volumetric  compaction.  At  high  pressures,  the  softening 

p 

term  goes  to  zero,  i.e.,  as  -*•  Wg  and  exp  (-Dlg  J£)  0  ,  the 

apparent  modulus  K  approaches  the  elastic  bulk  modulus,  K  . 

S  8 

Note  that  if  a  sample  is  first  isotropically  consolidated  (from 
point  1  to  point  2  in  Figure  5.1),  then  unloaded  (point  2  to  point  3), 
and  then  reloaded  (point  3  to  point  2),  the  model  dictates  that  the 
unloading-reloading  behavior  is  purely  elastic. 

5.2.2  Undrained  Condition  | 

During  an  undrained  hydrostatic  loading,  the  effective  stresses  and 
the  pore  pressure  are  generally  not  zero.  The  stress-strain  relations 
for  this  drainage  condition  can  be  obtained  in  a  manner  similar  to  that 
used  to  derive  Equations  5.11  through  5.13,  except  that  these  new 
expressions  will  involve  the  total  stresses  (Equations  5-5  and  5.6) 
instead  of  the  effective  stresses;  i.e.,  the  model  coefficients  will 
have  the  subscript  m  (mixture)  instead  of  the  subscript  s  (skeleton) 
(Table  fc.l);  thus: 
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5.2.3  Computation  of  Effective 
Stress  and  Pore  Pressure 

The  effective  stresses  and  the  pore  pressures  generated  during 
undrained  hydrostatic  loading  tests  (in  which  the  applied  total  stresses 
are  known)  can  be  computed  using  the  assumption  that  the  volumetric 
strains  from  Equations  5.13  and  5-17  are  equal  (Procedure  1  of  Chap¬ 
ter  U ) ;  thus : 
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from  which  can  be  obtained  as  a  function  of  using  a  bit  of 
algebra.  The  effective  stresses  and  the  pore  pressure  u  then  become 
(see  Equation  2.4) 


o' 

z 


(5.19) 


u 


(J1  -  Ji> 


(5.20) 


Note  that  when  the  material  is  fully  saturated  (i.e.,  a  two-phase 
system)  and  the  water  is  assumed  to  be  incompressible,  the  right  side  of 
Equation  5.17  or  5.18  becomes  zero  (i.e.,  is  independent  of  J  ), 
and  Equation  5.18  can  be  satisfied  if  and  only  if  is  equal  to  zero. 
This  means  that  all  of  the  applied  load  is  carried  by  the  water. 

5.3  SHEAR  PHASE 

During  the  shear  phase  of  a  conventional  TX  test,  the  cell  pressure 
is  held  constant  while  the  axial  stress  is  changed;  i.e.. 


constant 


(5.21) 


and 


do  -  0  (5.22) 
r 

where  Pc  is  the  total  confining  pressure  at  the  end  of  the  hydrostatic 
compression  phase.  If  the  hydrostatic  compression  phase  preceding  shear 
was  drained,  the  confining  pressure  Pc  is  also  the  effective  confining 
pressure  P^  . 

5.3.1  Drained  Condition 

During  a  drained  shear  test,  the  effective  and  total  stresses  are 
equal  (i.e.,  the  effective  stress  path  is  known  and  is  identical  with 
the  total  stress  path). 


f*  r'SPOnS'  °f  ““  8011  8ltelrt°n  *•  using  Equation  A.  29  of 
ppendix  A  and  the  response  function,  and  material  parameters  listed 
under  Set  1  of  Table  h.l  (for  the  drained  condition).  Thus: 
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for  the  drained  condition  (Table  1»,1). 

In  vie.  of  Equations  A.13  and  A.22  of  Appendix  A,  the  deviatoric 
components  of  the  total  strain  increment  tensor  are 
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According  to  Equation.  A.22  and  3.3  or  3.1,  the  value  of  d>  can 
be  written  as  s 
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where 
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dig  *  ~  de  deA  *  second  invariant  of  the  plastic  strain 
J  “  increment  deviation  tensor 


Substitution  of  Equation  5.25  intc  Equation  5.2h  leads  to 
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For  TX  test  conditions,  Equation  5.26  can  be  written  as 
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where  G  is  the  apparent  shear  modulus  of  the  material  under  drained 
s 

triaxial  loading  conditions.  The  second  term  on  the  right  3ide  of 
Equation  5.28  produces  a  softening  of  the  apparent  loading  shear  modulus 
due  to  plastic  flow. 

Within  the  yield  surface,  the  plastic  strain  increments  are  zero 
and  the  apparent  shear  modulus  equals  the  elastic  shear  modulus  G^  . 

The  volumetric  strain  increment  can  be  obtained  by  multiplying  both 
sides  of  Equation  5.23  by  the  Kronecker  delta,  6^  ,  and  then  using 
Equation  5*25;  thus: 
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Typical  (qualitative)  results  predicted  by  the  new  model  for  two 
types  of  drained  TX  shear  tests  are  shown  in  Figure  5.2. 

5.3.2  Undrained  Conditions 

During  an  undrained  shear  test  (following  isotropic  consolidation) 
only  the  total  stress  path  and,  consequently,  the  toted  stresses  are 
known  for  the  skeleton- water-air  mixture.  The  material  response  for 
this  test  can  be  modeled  using  the  mixture  model  from  Table  U.l,  and 
equations  similar  to  Equations  5.23  through  5.29  can  be  developed  except 
that  they  will  involve  the  total  stresses  (Equations  5.5  and  5.6) 
instead  of  the  effective  stresses.  These  model  coefficients  will  have 
the  subscript  m  (mixture)  instead  of, the  subscript  s  (skeleton). 

The  resulting  total  strain  increments  are 
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and  the  deviatoric  and  volumetric  strain  increments  are 


ij  2C„  V 


m 


'H 


ra 


where 


m 


F  (J.  ,  VjI  ,  k  )  on  the  hardening  surface 
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The  isotropic  compression  phase  preceding  shear  could  either  be 
drained  or  undrained.  In  this  example,  however,  it  was  chosen  to 
be  drained. 
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Equation  A. 28  of  Appendix  A  defines  dX  ,  and  consequently  (vdl„)  , 

m  '  2'a 

using  the  material  parameters  for  the  undrained  condition  (Table  U.l). 
Equation  5.31  can  be  written  in  a  form  similar  to  that  of  Equation 

5.28: 
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If  it  is  assumed  that  strain  difference  increments  are  the  same  under 
both  drained  and  undrained  conditions  and  hence  that  the  effective 
stress  paths  for  both  conditions  are  identical,  then  Equations  5.28  and 
5.33  lead  to 
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Equation  5-31*  dictates  that  the  elastic  moduli  of  the  material  under 
both  drained  and  undrained  shear  test  conditions  will  be  the  same  if  the 
material  undergoes  the  same  plastic  flow  under  both  conditions;  i.e., 

lf  (^l  ■  (^11  •  «<•»  • 

5.3.3  Effective  Stress  and  Pore 
Pressure  Computations 

During  the  undrained  shear  test,  the  effective  stresses  and  the 
pore  pressure  at  the  end  of  each  applied  total  stress  loading  increment 
can  be  computed  through  the  assumption  that  the  total  volumetric  strain 
increments  obtained  from  Equations  5.29  and  5 .32  are  equal  (i.e.. 
Procedure  1_ of  Chapter  U).  This  assumption  entails  the  added  assumption 
that  tne  pore  fluid  has  no  effect  on  the  shear  behavior  of  the  material. 
The  strain  compatibility  constraint  gives 
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from  which  can  be  obtained  as  a  function  of  .  The  pore  pres¬ 
sure  and  the  effective  stresses  then  become  (see  Equation  2.4) 
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Typical  (qualitative)  results  predicted  by  the  new  model  for  undrained 
TX  shear  tests  are  shown  in  Figure  5.3.  Figure  5.3  also  depicts  quali¬ 
tatively  the  effects  of  the  parameter  R  on  the  stress-strain  and  pore 
pressure  responses  during  a  conventional  undrained  shear  test. 

When  the  material  is  fully  saturated  (i.e.,  a  two-phase  system) 
and  the  water  is  assumed  to  be  incompressible,  the  right-hand  side  of 
Equation  5*35  becomes  zero;  i.e.,  is  independent  of  .  This 
means  that  the  effective  stress  path  is  independent  of  the  total  stress 
path  applied  to  the  material.  This  behavior  is  predicted  by  the  model 
reported  in  Reference  15> 


total  stress 
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CHAPTER  6 


COMPARISONS  OP  LABORATORY  TEST  DATA  WITH  MODEL  PREDICTIONS 

6.1  EXPERIMENTAL  PROGRAM 

The  material  used  in  the  experimental  program  was  Campbell  Svamp 
sand.  A  typical  gradation  curve  for  ihis  sand  is  shown  in  Figure  6.1. 

The  Campbell  Swamp  sand  is  classified  as  a  tan  uniform  sand  (SP)  and 
consists  of  subrounded  to  subangular  particles.  The  experimental 
program  consisted  of  a  total  of  52  TX  tests.  Eighteen  of  these  were 
conducted  on  15  percent  saturated  specimens  (designated  as  moist  tests), 
twenty-three  tests  were  conducted  on  98  percent  saturated  specimens 
(designated  as  vet  tests),  and  eleven  tests  were  conducted  on  100  per¬ 
cent  saturated  specimens  (designated  as  saturated  tests).  The  degree  of 
saturation  in  the  moist  and  vet  specimens  vas  achieved  by  mixing  water 
and  sand  before  building  the  specimens.  The  100  percent  saturation  of 
the  saturated  specimens,  however,  vas  achieved  by  back-pressure  satura¬ 
tion  (BPS)  techniques.  Based  on  their  degree  of  saturation,  the  52  test 
specimens  vere  naturally  divided  into  three  groups. 

The  first  group  of  specimens  (moist)  had  an  average  dry  density1  of 

3  3 

l.«?  g/cm  and  an  average  vet  density  of  1.^9  g/cm  .  The  folloving  tests 

vere  conducted  on  this  group:  load/unload  drained  isotropic  compression 

(consolidation)  tests,  load/unload  drained  uniaxial  strain  tests,  and 

load/unload  consolldated-drained  triaxial  shear  tests. 

The  second  group  of  specimens  (vet)  bad  an  average  dry  density  of 
1.62  g/cm}  and  an  average  vet  density  of  1.99  g/cm^.  This  group  vas 
tested  as  follows:  load/unload  drained  and  un drained  isotropic  compres¬ 
sion  tests,  load/unload  drained  and  undrained  uniaxial  strain  tests,  and 
load/unload  consolidated-drained  and  consolldated-undr mined  triaxial 
shear  tests  vith  pore  pressure  measurements  taken  on  all  the  undralned 
tests. 

1  A  table  of  factors  for  converting  metric  (Si)  units  of  measurement  to 
U.  S.  customary  units  and  U.  S.  customary  units  to  metric  (SI)  units 
is  found  on  page  7. 
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The  third  group  of  specimens  (saturated)  had  an  average  dry  density 
of  1.62  g/cm3  and  an  average  vet  density  of  2.01  g/em3.  This  group  was 
tested  as  follows:  load/unload  drained  isotropic  compression  teats  and 
load/unload  consolidated-drained,  consolidated-undrained,  and 
unconsolidated-undrained  triaxial  shear  tests,  with  pore  pressure  mea¬ 
surements  taken  on  the  undrained  tests. 

The  complete  experimental  program  is  summarized  in  the  following 
tabulation: 


Experimental  Program 

Group  1 
(Moist) 

Group  2 
(Wet) 

Group  3 
(Saturated) 

3 

Dry  density,  g/cm 

l.*»j 

1.62 

1.62 

Wet  density,  g/cm 

1.1*9 

1.99 

2.01 

Degree  of  saturation,  % 

15.0 

98.0 

100.0 

Number  of  tests  conducted 

18 

23 

11 

rIC 

IC 

IC 

DUX 

UIC 

CDTX 

CDTX 

DUX 

CUTX 

•type  of  tests  conducted11 

< 

UUX 

UUTX 

CDTX 

CUTX 

a 


C1VTX  -  Consolidated-drained  triaxial  shear  test 
CUTX  -  Consolidated-undrained  triaxial  shear  test 
DUX  -  Drained  uniaxial  strain  test 
IC  -  Isotropic  compression  (consolidation)  test;  drained 
UIC  -  Undrained  isotropic  compression  test 
UUTX  -  Unconsolidated-undrained  triaxial  shear  test 
UUX  -  Undrained  uniaxial  strain  test 


6.2  MATERIAL  CONSTANTS 

As  indicated  in  Chapter  *«,  a  maximum  of  *»2  material  constants  are 
associated  with  the  proposed  multiphase  constitutive  model  (see  Table 
*».l),  which  must  be  determined  experimentally  by  fitting  the  model  to 
laboratory  test  results.  Twenty-one  of  the  material  constants  are 
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associated  with  the  drained  behavior,  and  twenty-one  are  associated  with 
the  undrained  behavior.  The  actual  number  of  constants  needed  to 
simulate  the  behavior  of  a  given  material,  however,  depends  on  its 
degree  of  saturation.  If  the  material  is  partially  saturated  (i.e.,  the 
wet  group  of  Campbell  Swamp  sand  tests),  all  42  constants  are  generally 
needed.  If  the  material  is  fully  saturated,  fewer  constants  are  needed. 
For  example,  for  the  Group  3  Campbell  Swamp  sand,  only  24  material 
constants  were  needed  (Tables  6.1  and  6.2);  21  of  these  are  associated 
with  the  drained  behavior  (subscript  s  ,  Table  U . 1 ) ,  while  3  are  asso¬ 
ciated  with  the  undrained  isotropic  compression  behavior  (subscript  m). 
Finally,  if  the  material  is  nearly  dry  (i.e.,  the  water  content  is  negli¬ 
gible,  such  as  occurs  in  the  moist  group  of  Campbell  Swamp  sand  tests), 
the  pore  air  pressures  should  be  small;  hence,  only  21  material  constants 
are  needed.  These  are  the  material  constants  listed  in  Table  4.1  with 
the  subscript  s  . 

The  numerical  values  of  the  material  constants  for  the  moist 
(Group  l),  the  wet  (Group  2),  and  the  saturated  (Group  3)  Campbell  Swamp 
sand  tests  are  given  in  Table  6.1  for  the  drained  behavior  (subscript  s) 
and  Table  6.2  for  the  undrained  behavior  (subscript  m).  The  deriva¬ 
tion  of  the  numerical  values  of  these  constants  is  not  discussed  herein. 
The  purpose  of  this  chapter  is  simply  to  elucidate  the  capabilities  of 
the  new  model. 

6.3  COMPARISON  OF  TEST  RESULTS  WITH  MODEL  BEHAVIOR 

Figures  6.2  through  6.4  compare  actual  teBt  results  with  the  model 
fits  for  the  moist  (Group  l)  material.^-  Figure  6.2  compares  effective 
mean  normal  stress  versus  volumetric  strain  for  the  isotropic  compres¬ 
sion  test.  Effective  stress  paths  for  uniaxial  strain  and  the  effec¬ 
tive  failure  envelope  are  compared  in  Figure  6.3.  Figure  6.4  compares 


Because  the  initial  degree  of  saturation  of  this  group  of  testa  was 
only  15  percent,  no  undrained  tests  were  conducted  on  this  material. 
However,  it  is  anticipated  that  the  undrained  behavior  of  these  tests 
could  also  be  represented  by  Figures  6.2  through  6.4. 
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consolidated  drained  triaxial  shear  test  response  for  various  confining 
pressure  levels. 

Test  results  are  compared  with  model  behavior  for  the  vet  (Group  2) 

2 

material  in  Figures  6.5  through  6.8.  Figure  6.5  compares  effective  and 
total  mean  normal  stress  versus  volumetric  strain.  Effective  and  total 
stress  paths  in  drained  and  undrained  states  of  uniaxial  strain  are  com- 
ared  in  Figure  6.6.  Figure  6.7  compares  consolidated  drained  triaxial 
shear  test  response.  The  results  of  consolidated  undrained  triaxial 
shear  tests  are  compared  with  the  model's  behavior  in  Figure  6.8. 

Figures  6.9  through  6.13  compare  test  results  versus  model  behavior 
for  the  saturated  (Group  3)  material.  The  comparison  plots  for  iso¬ 
tropic  consolidation  are  shown  in  Figure  6.9-  Figure  6.10  compares  con¬ 
solidated  drained  triaxial  shear  response.  Figures  6.11  and  6.12  com¬ 
pare  consolidated  undrained  triaxial  shear  response.  Principal  stress 
difference  versus  axial  strain  and  pore  pressure  versus  axial  strain  for 
unconsolidated  undrained  triaxial  shear  are  compared  in  Figure  6.13. 

Figures  6.2  through  6.13  clearly  show  that  the  proposed  constitu¬ 
tive  model  qualitatively  simulates  the  various  stress-strain  and  pore 
pressure  responses  of  Campbell  Swamp  sand  for  the  moist  (Group  l),  the 
wet  (Group  2),  and  the  saturated  (Group  3)  conditions.  It  should  be 
pointed  out  that  pore  pressure  measurements  are  not  needed  for  fitting 
the  new  model;  they  are  useful,  however,  for  verifying  the  resultant 
model  fits. 


2 

The  drained  and  undralned  tests  conducted  on  the  vet  (Group  2)  sand 
produced  identical  results;  i.e.,  no  pore  pressures  were  measured  in 
the  undrained  tests.  It  is  anticipated  that  at  stress  levels  above 
those  employed  in  the  current  experimental  program  (20  bars)  the  wet 
sand  would  approach  and  eventually  reach  full  saturation  during 
undrained  tests  and  consequently  the  undrained  test  results  would 
-  diverge  from  the  drained  test  results. 

The  bulk  modulus  of  this  material  under  undralned  isotropic  compression 
was  assumed  to  be  constant  and  equal  to  1*5  kbars  (see  Table  6.2). 

This  value  was  calculated  from  mixture  theory;  no  undrained  isotropic 
compression  test  was  actually  conducted  on  this  particular  material. 
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Table  6.1.  Numerical  values  of  material  constants  for 
drained  behavior  of  Campbell  Swamp  sand. 


Material  Constants 

Group  1 

Group  2 

Group  3 

Name 

Notation,  Units 

(Moist) 

(Wet) 

( Saturated ) 

Dry  Density 

Ps»  g/cm3 

1.1*3 

1.62 

1.62 

A  ,  kbars 
s 

105.0002886 

120.0001 

0.288675 

Failure  Envelope  4 

B  ,  (kbars)  1 

0.0023 

0.0023 

0.96 

Parameter 

s 

C  ,  kbars 

l  3 

105.0 

120.0 

0.288669 

f  Ris*  " 

3.1 

3.1 

3.1 

R,  ,  — 

0.0 

0.0 

0.0 

Is  -1 

R2g,  (kbars) 

0.0 

0.0 

0.0 

R,  ,  — 

-1.0 

0.0 

0.0 

33  „p 

R^s ,  (kbars) 

R_  ,  kbars 

3000.0 

0.0 

0.0 

Hardening  Surface  4 

0.011 

0.0 

0.0 

Parameters 

5s 

W  ,  — 

3 

0.305 

0.325 

0.325 

D  ,  (kbars)”1 
s 

1.825 

0.23 

0.23 

W,  ,  (kbars)”2 
13  -1 

55.0 

39.0 

39.0 

D.  ,  (kbars) 

120.0 

130.0 

130.0 

^  Is 

r  K.  ,  kbars 
is 

0.1* 

0.5 

0.5 

Bulk  Modulus  4 

“is*  ~  -1 

K  ,  (kbars) 

2s 

0.92 

0.75 

0.75 

Parameters 

55.0 

60.0 

60.0 

f  Ci.  ,  kbars 
is 

0.2 

0.3l* 

0.3l* 

°ls’  ""  -1 

0.95 

0.75 

0.75 

Shear  Modulus  4 

G2g,  (kbars)” 

G_  ,  kbars 

j3 

200.0 

21*0.0 

21*0.0 

Parameters 

2.0 

1.0 

1.0 

*.  G,  ,  — — 
l*S 

200.0 

10.0 

10.0 
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Table  6.2.  Numerical  values  of  material  constants  for 
undrained  behavior  of  Campbell  Swamp  sand. 


Material  Constants  Group  1  Group  2  Group  3 

Name _  Notation.  Units  (Moist)  (Wet)  ( Saturated ) 


Wet  Density 


p  ,  g/cm3 
m 


Failure  Envelope 
Parameters 


Hardening  Surface 
Parameters 


kbars 

NA 

120.0001 

0.2886 

(kbars)  1 

NA 

0.0023 

0.96a 

kbars 

NA 

120.0 

0.2886 

— 

NA 

3.1 

NAb 

— 

NA 

0.0 

NAb 

(kbars )-1 

NA 

0.0 

NAb 

— 

NA 

0.0 

NAb 

(kbars )”2 

NA 

0.0 

NAb 

kbars 

NA 

0.0 

NAb 

— 

NA 

0.016 

NAb 

(kbars)”1 

NA 

1.0 

NAb 

(kbars)"2 

NA 

30.0 

NAb 

(kbars)  1 

NA 

95.0 

NAb 

Bulk  Modulus 
Parameters 


Shear  Modulus 
Parameters 


K,  ,  kbars 
ini 


K2m>  (kbars)" 


G.  ,  kbars 
ira 


C2m.  (kbars)" 
G^,  kbars 


0.5 

1*5.0 

0.99 

0.0 

1.30 

0.0 

60.0 

0.3l*  H 

0.3U 

0.75a 

0.75 

2l*0.0a 

!  1*0.0 

1.0° 

10.0 

10.0a 

The  failure  envelope  parameters  and  the  shear  modulus  parameters  for 
this  material  are  the  same  as  those  in  Table  6.1. 

For  this  group  the  hardening  surface  was  not  utilized  in  the 
calculations. 
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Figure  6.2 
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Drained  isotropic  compression  response, 
ments  versus  model  behavior.  Group  1. 


laboratory  measure- 
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Figure  6. U  Principal  stress  difference  versus  axial  strain  response 
from  consolidated  drained  triaxial  shear  tests,  labora¬ 
tory  measurements  versus  model  behavior.  Group  1. 
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Figure  6.5  Drained  and  undralned  Inotropl.-  ompreaslon 
response,  laboratory  measurements  versus 
nodel  behavior.  Group 
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EFFECTIVE  AND  TOTAL  MEAN  NORMAL  STRESS  P  AND  P,  BARS 

6.6  Drained  and  undrained  stress  path  for  uniaxial  (UX)  strain,  laboratory  measurements 
versus  model  behavior.  Group  2. 


principal  stress 


i  a 
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AXIAL  STRAIN  «r  ,  PERCENT 

Figure  6.8  Principal  stress  difference  and  excess  pore 
pressure  versus  axial  strain  response  from 
consolidated  undrained  triaxial  shear  tests, 
laboratory  measurements  versus  model 
behavior.  Group 
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Figure  6.9  Drained  iaotropic  coepression  response,  laboratory 
measurements  versus  model  behavior.  Group  3. 
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Figure  to.  10  Principal  stress  difference  versus  axial 
strain  response  Trom  consol  1 dated 
drained  triaxlal  shear  teats,  laboratory 
measurements  versus  model  behavior, 

Croup  3. 
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Figure  6.12  Principal  stress  difference  and  excess 
pore  pressure  versus  axial  strain 
response  from  consolidated  undrained 
triaxial  shear  tests,  laboratory 
measurements  versus  model  behavior. 
Group  3. 
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Figure  6,13  Principal  stress  difference  and  excess  pore 
pressure  versus  axial  strain  response  from 
unconsolidated  undrained  triaxial  shear 
tests,  laboratory  measurements  versus  model 
behavior,  Group  3. 


CHAPTER  7 


APPLICATION  OF  THE  MODEL  FOR  SOLVING  GROUND  SHOCK  PROBLEMS 

7.1  INTRODUCTION 

The  multiphase  constitutive  model  described  in  Chapter  was 
incorporated  into  the  LAYER  code  to  provide  an  effective  stress  ground 
motion  calculation  capability  for  wet  sites.  To  demonstrate  this 
capability,  the  code  was  used  to  perform  a  2D  axisymmetric  calculation 
using  the  saturated  material  properties  and  model  fits  presented  in 
Chapter  6. 

7.2  COMPUTATIONAL  DETAILS 

The  basic  geometry  for  the  calculation  is  shown  schematically  in 
Figure  7.1.  The  soil  profile  consisted  of  a  half-space  of  the  saturated 
(Group  3)  material.  As  shown  in  Figure  7.1,  the  domain  of  the  calcula¬ 
tion  extended  in  the  radial  direction  for  a  distance  of  1*1.8  metres  from 
ground  zero  (G2)  and  vertically  to  a  depth  of  1*1.8  metres;  i.e.,  the 
region  of  the  calculation  was  square.  A  0.1* -metre- square,  finite- 
difference  grid  size  and  a  time  step  of  0.15  ms  were  used  to  carry  the 
2D  finite-difference  calculation  to  a  real  time  of  18  ms. 

The  far  right  and  bottom  boundary  conditions  were  of  the  transmit¬ 
ting  type  (Reference  20).  The  left  boundary  was,  of  course,  treated  as 
an  axisymmetric  boundary;  i.e.,  the  radial  velocity  and  the  shear  strain 
(or  stress)  were  both  required  to  vanish.  The  surface  boundary  was 
subjected  to  a  radially  expanding  airblatt  loading,  starting  at  GZ  at 
zero  time  (Figure  7.2),  The  airblast  shock  velocity  VS  was  constant 
and  equal  to  3  metres  per  millisecond.  The  peak  overpressure  Pmny 
was  also  constant  and  equal  to  100  bars.  A  constant  rise  time  t  and 

a  constant  positive  phase  duration  d  of  3  ms  and  9  ms,  respectively, 
were  used  to  characterize  the  shape  of  the  airslap  pulse. 

The  far  right  and  bottom  transmitting  boundary  conditions  were 
applied  at  velocity  points  while  the  left  and  surface  boundary  condi¬ 
tions  were  applied  at  stress  points. 
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7.3  RESULTS 

In  this  section,  selected  code  stress-  and  motion-time  histories 
and  stress  path  plots  are  presented  at  two  intermediate  ground  ranges: 

R  *  3.6  metres  and  15.6  metres.  These  ranges  pertain  to  the  stress 
output;  for  motion,  the  corresponding  ranges  are  R  *  3.8  metres  and 
15.8  metres.  The  depths  are  1.6,  3.6,  and  7.6  metres  for  stresses  and 
1.8,  3.8,  and  7.8  metres  for  motions. 

Figures  7.3  through  7.5  present  total  vertical  stress,  total  radial 
stress,  total  mean  normal  stress,  second  invariant  of  the  stress  devia¬ 
tion  tensor,  pore  pressure,  and  effective  mean  normal  stress  for  the 
3.6-metre  range  and  the  1.6-,  3.6-,  and  7.8-metre  depths.  Figures  7.6 
through  7.8  show  the  corresponding  plots  for  the  15.6-metre  range.  It 
is  clear  from  these  figures  that  the  effective  stresses  are  very  small 
in  comparison  with  the  total  stresses  and  that  most  of  the  stress  is 
carried  by  the  pore  fluid.  This  was  expected  because  the  material  is 
fully  saturated. 

Figures  7.9  and  7-10  present  the  vertical  and  radial  particle 
velocity  and  displacement  wave  forms,  respectively,  for  the  3.8-metre 
range  at  the  1.8-,  3.8-,  and  7.8-metre  depths.  Corresponding  motion 
wave  forms  for  the  15.8-metre  range  are  shown  in  Figures  7-11  and  7.12. 
The  radial  velocities  and  displacements  are  smaller  than  the  vertical 
motions  because  the  problem  is  superseismic. 

Stress  path  plots  for  the  3.6-metre  range  and  1.6-,  3.6-,  and  7.6- 
metre  depths  are  presented  in  Figure  7.13.  Figure  7.1^  depicts  the 
corresponding  plots  for  the  15.6-metre  range. 
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Figure  7«1  Axi symmetric  problem  geometry. 
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Figure  7.6  Stress-tine  histories;  range  15.6  metres,  depth  1.6  metres 


P',  BARS 

P  ,  BARS 

a.  DEPTH  1.6  METRES 

P,  BARS 

P  ,  BARS 

b.  DEPTH  3.6  METRES 

•  0 

P  ,  BARS 

C.  DEPTH  7.6  METRES 


Figure  7.13  Comparison  of  total  and  effective  stress  paths;  range 
3.6  metres. 
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CHAPTER  8 


SUMMARY  AND  RECOMMENDATIONS 


8.1  SUMMARY 

A  three-dimensional,  elastic-plastic,  work-hardening  constitutive 
model  for  simulating  isotropic  three-phase  earth  materials  has  been 
developed.  Within  the  elastic  range  the  model  contains  two  response 
functions:  a  nonlinear  elastic  bulk  modulus  and  a  nonlinear  elastic 
shear  modulus.  In  the  plastic  range  it  contains  three  potential  func¬ 
tions:  a  failure  envelope,  a  hardening  yield  surface,  and  a  hardening 
function.  The  numerical  values  of  the  coefficients  in  these  response 
and  potential  functions  can  be  determined  experimentally  using  essen¬ 
tially  conventional  soil  testing  techniques.  The  major  feature  of  the 
new  constitutive  model  is  its  capability  of  simulating  the  total  and 
effective  stress-stain  and  pore  pressure  responses  of  saturated  and 
partially  saturated  soils. 

The  behavior  of  the  model  under  drained  and  undrained  triaxial  test 
conditions  has  been  examined,  and  the  method  for  calculating  the  effec¬ 
tive  and  total  stresses  and  the  pore  pressure  of  a  deforming  multiphase 
soil  subjected  to  given  total  stress  or  strain  paths  has  been  explained. 

The  new  model  has  been  fit  to  a  set  of  drained  and  undrained 
mechanical  property  test  data  for  moist,  wet,  and  saturated  Campbell 
Swamp  sand;  and  the  saturated  sand  model  has  be*n  used  in  a  demonstra¬ 
tion  LAYER  code  calculation  to  show  its  potential  for  performing  effec¬ 
tive  stress  ground  shock  analyses.  The  demonstration  calculation  was 
precisely  what  its  name  implied,  a  demonstration.  The  problem  that  was 
solved  was  trivial;  the  effort  to  solve  it,  however,  was  not. 

8.2  RECOMMENDATIONS 

It  is  recommended  that  the  new  model  be  verified  through  calcula¬ 
tions  against  a  wide  variety  of  field  events  executed  at  wet  sites. 

These  efforts  should  also  include  calculations  with  current  single-phase 
(total  stress)  constitutive  models  so  that  the  benefits  of  effective 
stress  ground  shock  analyses  can  be  assessed. 
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APPENDIX  A 

FUNDAMENTAL  BASIS  OF  ELASTIC-PLASTIC  CONSTITUTIVE  MODELS 

A.l  BASIC  CONCEPTS  FROM  CONTINUUM  MECHANICS 

From  a  microscopic  point  of  view,  physical  bodies  are  composed  of 
discrete  molecules  interconnected  by  internal  forces  of  mutual 
attraction  and  repulsion.  The  concept  of  stress  within  a  body  requires 
that  boundary  distances  and/or  loaded  area  be  large  in  comparison  with 
distances  between  molecules  and/or  the  size  of  the  individual  molecule. 
This,  in  effect,  transforms  a  body  composed  of  discrete  molecules  into  a 
statistically  macroscopic  equivalent  amenable  to  mathematical  analysis. 
Since  most  engineering  problems  deal  with  macroscopic  phenomena  and 
involve  very  large  boundary  distances  and  loaded  areas  compared  with 
individual  molecules,  it  appears  reasonable  and  convenient  to  invoke  the 
mechanics  of  continua  as  the  basis  for  analytical  consideration  of  these 
problems . 

The  theory  of  continuous  meuia  is  built  upon  two  strong  founda¬ 
tions:  basic  balance,  and  conservation  laws  and  a  constitutive  theory. 

The  basic  balance  and  conservative  lavs  of  any  continuum  are 

1.  Conservation  of  mass 

2.  Conservation  of  energy 

3.  Balance  of  linear  momentum 
b.  Balance  of  angular  momentum 

5.  Inadmissibility  of  decreasing  entropy 

When  thermal  effects  are  neglected  and  symmetry  of  the  stress  tensor  is 
assumed,  these  basic  axioms  of  continuum  mechanics  lead  to  the 
following  continuity  equation:1 

foV.i  ■ 0  (A*1) 

1  Indices  take  on  a  value  of  1,  2,  or  3.  A  repeated  index  is  to  be  summed 
over  its  range.  A  comma  between  subscripts  represents  a  derivative. 
Quantities  are  referred  to  rectangular  Cartesian  coordinates  . 


and  the  equations  of  motion: 

°lj,  j  *  Fi  -  pai  =  0  (A'2) 

where 

p  =  mass  density 
t  =  time 

v =  components  of  velocity  vector 
°ij  =  °Ji  *  sJr™ietr^cal  stress  tensor 
=  components  of  body  force 
a^^  =  components  of  acceleration  vector 

Equations  A.l  and  A. 2  are  called  field  equations.  They  constitute 
four  equations  that  involve  ten  unknown  functions  of  time  and  space. 
Therefore,  the  system  resulting  from  Equations  A.l  and  A. 2  is  indetermi¬ 
nate.  These  unknown  functions  sure  the  mass  density  p  ,  the  three 
velocity  components  v^  ,  and  the  six  independent  stress  components 
Ojj  .  The  body  force  components  ^  are  known  quantities,  and  the 
acceleration  components  ai  are  expressible  in  terms  of  the  velocity 
components  v^  .  To  overcome  the  indeterminacy  and  make  the  system 
complete,  six  additional  expressions  relating  the  ten  unknown  variables 
are  required.  In  continuum  mechanics  such  relations  are  stated  by 
constitutive  equations  (or  material  models),  which  relate  stresses  to 
deformation  and  history  of  deformation.  The  difference  between  consti¬ 
tutive  equations  and  field  equations  (Equations  A.l  and  A. 2)  is  that  the 
latter  contains  both  space  coordinates  and  time  and  is  applicable  to  all 
materials,  whereas  the  former  is  independent  of  space  coordinates  (for 
homogeneous  materials)  and  represents  the  intrinsic  response  of  a 
particular  material  or  class  of  materials  and,  as  such,  is  a  mathemati¬ 
cal  idealization  of  the  mechanical  behavior  of  real  materials. 

The  general  form  of  a  constitutive  equation  may  be  expressed  by  the 
following  functional  form: 
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HlJ(Dmn  ■  °w  •  Ers  *  °kl  *  p)  -  0 


(A. 3) 


where  cr^  is  the  components  of  stress  tensor  and  the  deformation-rate 

tensor  D  and  the  spin  tensor  fl  are  related  to  the  components  of 
mn  r  qp 

the  velocity  vector,  v^  ,  by 


D  *  (v  +  v  ) 
mn  2  m,n  n,m 


and 


(A. It) 


q  *  —  (v  -  v  ) 
OP  2  q.p  p,q 


and  th:  infinitesimal  strain  tensor  cr8  is  related  to  the  components 
of  the  displacement  vector  by 


c 


rs 


u 


s.r 


) 


(A. 5) 


Equations  A.l  through  A. 3  constitute  ten  equations  involving  ten  unknown 
variables.  These  equations  will  lead,  in  conjunction  with  the  kinematic 
relations  given  by  Equations  A. It  and  A. 5  and  the  appropriate  boundary 
conditions,  to  a  complete  description  for  the  solution  of  a  boundary- 
value  problem. 

In  general,  materials  having  the  same  mass  and  geometry  respond 
differently  when  subjected  to  identical  external  effects.  Therefore,  a 
variety  of  constitutive  theories  have  emerged,  each  of  which  describes  a 
limited  number  of  physical  phenomena  decided  on  at  the  outset  for  a 
given  material.  In  the  following  sections  the  constitutive  theory  for 
an  elastic-plastic  isotropic  material  is  presented. 

A. 2  GENERAL  DESCRIPTION  OP  ELASTIC-PLASTIC 

CONSTITUTIVE  MODELS 

The  basic  premise  of  elastic-plastic  constitutive  models  is  the 
assumption  that  certain  materials  are  capable  of  undergoing  small 
plastic  (permanent)  as  well  as  elastic  (recoverable)  strains  at  each 
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The  behavior  of  the  material  in  the  plastic  range  can  be  described  with¬ 
in  the  framework  of  the  generalised  incremental  theory  of  plasticity. 

The  mathematical  basis  of  the  theory  was  established  by  Drucker  (Refer- 
ence  21),  who  introduced  the  concept  of  material  stability,  which  has 
the  following  implications: 

1.  Yield  surface  (leading  function)  should  be  convex  in  stress 
space. 

2.  Yield  surface  and  plastic  potential  should  coincide  (which 
results  in  an  "associated"  flow  rule). 


o 

References  cited  in  this  appendix  are  included  in  the  References 
at  the  end  of  the  main  text. 


9* 


3.  Work  ’'softening”  should  not  occur. 

These  three  conditions  can  be  summarized  mathematically  by  the  following 
inequality 

doij  dE«  ^ 0  (A-8) 

These  conditions  allow  considerable  flexibility  in  the  choice  of  the 
form  of  the  loading  function,  ,  for  the  model,  which  serves  as  both 
a  yield  surface  and  the  plastic  potential.  In  general,  the  yield 
surface  may  be  expressed  as 

i  (ot1  ,  tc)  -  0  (A. 9) 


and  for  isotropic  materials  the  yield  surface  may  be  expressed,  for 
example,  as 


i  (J,  ,  ,  <)  -  0 

A  O 


(A. 10) 


where 


J,  ■  o  ■  first  invariant  of  the  stress  tensor 
x  nn 

J2  ■  |  S^j  8  ■  second  invariant  of  the  stress  deviation  tensor 

-  vJj/3)  *  stress  deviation  tensor 


■  Kronecker  delta 


(1  i  ■  J 

(0  i  +  J 


<  ■  a  hardening  parameter 


The  hardening  parameter,  >c  ,  generally  can  be  taken  to  be  a  function  of 

P 

the  plastic  strain  tensor,  .  The  yield  surface  of  Equation  A. 9  or 
A. 10  may  expand  or  contract  as  k  increases  or  decreases,  respectively 
(Figure  A.l). 
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Condition*  1  through  3  above,  taken  it.  conjunction  with  Equation 
A. 9  or  A. 10,  result  in  the  following  plastic  flow  rule  for  isotropic 

materials s 


if  {  ■  0 
if  i  <  0 


(A. 11) 


where  dX  is  a  positive  scalar  factor  of  proportionality,  which  is  non- 
sero  only  when  plastic  deformations  occur  and  is  dependent  on  the 
particular  form  of  the  loading  function. 

A. 2.1  Elastic  Strain  Increment  Tensor 

For  isotropic  elastic  materials,  the  strain  increment  tensor 
(Equation  A. 7)  takes  the  following  form 


The  bulk  and  shear  moduli  can  be  functions  of  the  invariants  of  the 
stress  tensor.  Accordingly,  it  is  assumed  that  K  ■  K(J^  ,  ,  J^} 

and  Q  «  ,  J^)  where  i*  the  third  invariant  of  the 

stress  deviation  tensor.  Equation  A.l?  can  be  written  in  terms  of  the 
hydrostatic  and  deviatoric  components  of  the  strain  and  stress  Increment 
tensors;  i.e,, 


(A. 13) 
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where 


E 

de^  ■  increment  of  elastic  volumetric  strain 
E 

de^  ■  elastic  strain  deviation  increment  tensor 


In  order  not  to  generate  ener w  or  hysteresis  within  the  elastic  ranpe, 
the  elastic  behavior  of  the  model  must  be  path  independent.  The  mate¬ 
rial  should  then  possess  a  positive  definite  elastic  internal  energy 
function  V  which  is  independent  of  stress  path.  The  strain  energy 
function  can  be  written  as 


W  - 


0 


♦iW 


Ji  Ji  sij  1 

-  r - — —  dJ  ♦  / J — —  s  ds 

0  9K(J  ,  J2  ,  J  )  1  0  ZGUj  ,  J2  ,  J3)  1J  1J 


J,  d(j  )2 
„  /A  - i - 

0  18K(J  ,  J 


(A.lU) 


In  order  for  V  to  be  independent  of  stress  path,  the  integrals  in 
Equation  A.lU  have  to  depend  only  on  the  current  values  of  and 
J2  .  Therefore,  the  bulk  and  shear  moduli  have* to  be  expressed  as 


K  -  K(JX) 
G  -  0(J2) 


(A. IS) 
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k 


Further,  K  and  G  must,  always  be  positive.  Since  during  elastic 
deformation  the  hardening  parameter  (<)  in  constant,  the  bulk  and 
shear  moduli  can  also  be  expressed  as 


K  ■  K(J:  ,  k) 


g(j2  ,  <) 


(A.16) 


A. 2. 2  Plastic  Strain  Increment  Tensor 

The  plastic  strain  increment  tensor  is  given  by  Equation  A. 11  where 
the  loading  function  j  is  given  by  Equation  A. 9  or  A. 10.  The  harden¬ 
ing  parameter  in  Equation  A. 9  or  A. 10  could  be  taken  as  being  equal  to 

p 

the  plastic  volumetric  strain  i  thus 


P 

*  *  "kk 


(A. 17) 


The  use  of  Equation  A. 17  will  allow  the  yield  surface  to  expand  as  well 
as  to  contract  (Figure  A.l).  The  plastic  loading  criteria  for  the  func¬ 
tion  6  are  given  as 


>  0  for  loading 
■  0  for  neutral  loading 
l <  0  for  unloading 


(A. 18) 


Because  dc^  ■  0  during  unloading  or  neutral  loading,  as  well  as  for 
$<  0  ,  Equations  A. 12  and/or  A. 13  are  used  to  determine  the  purely 
elastic  strain  changes.  The  prescription  that  neutral  loading  produces 
no  plastic  strain  is  called  the  continuity  condition.  Its  satisfaction 
leads  to  coincidence  of  the  elastic  and  clastic  constitutive  laws  during 
neutral  loading  (References  21  and  22). 

Like  the  elastic  behavior,  the  plastic  stress-strain  relation  can 
be  expressed  in  terms  of  the  hydrostatic  and  deviatoric  components  of 
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‘'I?',  *h'  Ch“n  of  to  th.  right-hand 

side  of  Equation  A. 11  yields 


dcf  -dX  + 

J  V 1  «  sau 


id,  *  dX  ( 4  +  __i_  _  e  \ 

1J  h  1J  «E*5TS«J 


(A.19) 


Multiplying  both  sides  of  Equation  A.19  by  6 


gives 


de£.  *  3  dX 
kk  3J 


(A. 20) 


fTlitorl0  of  th.  plastic  .train  incr«.nt  tenaor  d.P 


can  be  written  as 


d«?,  ■  «<,  -  i  d 4  * 


(A. 21) 


Substitution  of  Equations  A.l9  and  A.20  into  Equation  A. 


21  yields 


def  -  -1L  s 


(A. 22) 


In  order  to  na.  Bluation  *.19  or  Equations  *.20  and  *.22,  the 

prooortionalit,  factor  dX  nust  he  detemined.  tti.  can  he  accon- 

P  Shed  in  th.  following  »nncr.  Prom  Equation.  *.l0  and  *.17  th.  total 
derivative  of  {  becomes 


dtf  «  «u,  +  -i - 2L 


3JX  ““1 


2-C  sC 


’«  “u  *rfdcLm° 


(A. 23) 


Hff 


Ia  view  of  Equations  A. 13  and  A. 20,  Equation  A. 23  becomes 


3K  de, 


G  de 


kk  3J, 


11 


-ii-s 

11 


♦  3  dX 


1L  JtL 

3J.  -  P 
1  8ckk 


(A. 24) 


Substituting  Equation  A. 6  into  Equation  A. 24  results  in 


SKfdc^  -  d^)  14-  ♦  ~Z  (d. 

1  rj2 


-  def .  S 
1J  »«:  iJ 


-3  dX 


Itff  (A-25> 

1  5ctt 


or 


3K  de, 


kk  3J 


l£_  +  -2 - 2i_ 


$2  8^2 


3iJ  dCiJ 


3K  dc; 


kk  3J. 


*ij  dCiJ 


-  3  dX 


a L_li_ 


(A. 26) 


3J 


i  3  c: 


kk 


p  p 

Substituting  the  values  of  de*  k  and  de^  from  Equations  A. 20  and 
A. 22,  respectively,  into  Equation  A. 26  yields 


3K  dc 


kk  3J. 


IL  +  _5__lL  s  de 
1  ^  iJ  iJ 


9K 


dX  (|f-)2  +  G  dX  (-^-) 


-  3  dX 


‘  3J 


aL  Ji_ 

3J1  3cP 
1  3ckk 


a/T 


(A. 27 ) 


Solving  for  dX  gives 
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dX 


3K  |4-  de, .  +  —  -ii_ 
1 


W,  “‘kk  T  ~  S, ,  de 


^  iJ  « 


9K 


♦  G  (•&  \  _  3  li _ i 

vi/  \>srj  >jiu 


SLJL. 

p 

kk 


(A. 28) 


A. 2. 3  Total  Strain  Increment  Tenser 

The  total  strain  increment  tensor  can  be  obtained  by  combining 
Equations  A.6,  A.12,  A.19,  and  A.28: 


dj.  as 

de .  »  6  +  — LI 

iJ  9K  iJ  2G 

l£_d, 

3J,  4ekk 


3K|^-de..  +-2--li-s  d 

/rjj/r  m  “ 


f~uy  .liLA L 


3J,  .  P 
1  3ckk 


(a. 29) 


/aj"6ii  +  — -^-s  \ 

\ 1  J  1Jj 


Similarly,  the  stress  increment  tensor  can  be  written 


as 


^ij  "  K  dcwu  ,  ♦  2G  de 


kk  “iJ 


'iJ 


f  3K  |4-  dc._  ♦  —  -24- 


3J,  kk 


_ S  de 

a/T  001  ““ 


J 


(A.30) 


f3K  |4-  fi.  ♦  g  -24 — LL 
1  ^ 


} 


Equation  A.29  or  Equation  *.30  1.  th.  «.n.ral  constltutlr.  squatlon 
r  elastic-plastic  Isotropic  aatsrlal.  To  uts  slthar  of  those 
aquations  It  1,  Mcecry  only  to  spwlfy  ths  functional  fonu  of  E 

.  aao  <  and.  of  course,  to  determine  erperinentelly  the  nu»rlcal’ 
▼mlues  of  the  coefficients  in  these  functions. 
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APPENDIX  B:  NOTATION 


Components  of  acceleration  vector 

A,  B,  and  C  Material  constants;  appear  in  the  equation  of  the 
failure  envelope 

^ijkl  Elastic  compliances  of  the  material 

d  Phase  duration;  appears  in  application  of  the  model 

dX  Positive  scalar  factor  of  proportionality;  appears 
in  the  flow  rule 

Dt  Dj  Material  constants;  appear  in  the  hardening  function 

D  Deformation-rate  tensor 

mn 

e^j  Strain  deviation  tensor 

g 

e^  j  Elastic  strain  deviation  tensor 

r 

Plastic  strain  deviation  tensor 

f  Failure  envelope 

4  Loading  function 

F  Strain-hardening  elliptic  cap 

Fj  Components  of  body  force 

G  Elastic  shear  modulus 

G.  Initial  elastic  shear  modulus 

k 

,  C^,  G^  Material  constants  in  shear  modulus  formulation 

l!.  Second  invariant  of  the  plastic  strain  deviation 
tensor 

First  invariant  of  the  stress  tensor 

Jj  First  invariant  of  the  effective  stress  tensor 

,  J.,  Second  invariant  of  the  effective  or  total  stress 
deviation  tensor*  respectively 

Third  invariant  of  the  str-rs  deviation  tensor 

K  Elastic  bulk  modulus 

Initial  elastic  bulk  modulus 
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*1**2 

L(k) 


■ 

P 

P' 

P 

• 

P 

c 

P' 

c 


r,0,« 

R 

RD 


t 


u 


u 


t 


V.. 


Material  constants  in  bulk  modulus  formulation 

Intersection  of  the  hardening  surface  with  the 
failure  envelope  f(J^ 

Subscript  used  to  indicate  "mixture" 

Total  aean  normal  stress 
Effective  mean  normal  stress 
Pore  air  pressure 

Total  confining  pressure  at  the  end  of  the 
hydrostatic  concession  phase  of  a  triaxial  test 

Effective  confining  pressure  at  the  end  of  the 
hydrostatic  compression  phase  of  a  triaxial  test 

Pore  water  pressure 

Cylindrical  coordinate  system 

Ratio  of  the  major  to  minor  axes  of  the  elliptic 
hardening  surface;  ground  ran;e 

Radial  displacement 

Material  constants  in  hardening  surface  formulation 

Subscript  used  to  indicate  "skeleton" 

Stress  deviation  tensor 

Time 

Rise  time 

Total  pore  pressu-e  repres  nting  the  combined  effect 
of  the  pore  air  pressure  and  the  pore  water  pressure 

Dlsplaceswnt  vector 

Components  of  the  velocity  vector 


U  Radial  particle  velocity 

V  Vertical  particle  velocity 

VD  Vertical  displacement 

VS  Alrblast  shock  velocity 

V  Maximus  plastic  volumetric  compression  that  the 
material  can  experience  under  hydrostatic  loading 

W  Elastic  internal  energy  function 


10* 


Material  constant 


Xj  Cartesian  coordinate  system 

X(k)  Intersection  of  the  hardening  surface  vlth  the 

Jx  axis 

Kronecker  delta 

AV/V  Volumetric  strain 

o 

c . .  Total  strain  tonsor 

E 

Cjj  Elastic  str&ln  tensor 

Total  volumetric  strain 

£ 

e. .  Elastic  volusetric  strain 

kk 

p 

Cjj  Plastic  strain  ensor 

p 

Plastic  volumetric  strain 

cr  Total  radial  strain 

c  Total  vertical  strain 

e0  Total  tangential  strain 

9  Mass  density 

Ojj  Total  stress  tensor 

Ojj  Effective  stress  tensor 

or  Total  radial  stress 

k  Hardening  parameter 

o'  Effective  radial  stress 

r 

os  Total  axial  stress 

o'  Effective  axial  stress 

£ 

0^  Total  tangential  stress 

o£  Effective  tangential  stress 
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Naval  Facilities  Engineering  Coaaumd 
Washington,  D.  C.  20390 
ATTN:  Code  09M22C  (Tech  Lib) 

Superintendent 
Naval  Postgraduate  School 
Monterey,  Calif.  93940 
ATTN:  Code  0142  (Library) 

Officer- in-Charge 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  Md.  20910 
ATTN:  Technical  Library 

Coaaander 

Naval  Surface  Weapons  Center 
Dahlgren,  Va.  22448 

ATTN:  Technical  Library  &  Information  Services  Branch 

Office  of  Naval  Research 
Arlington,  Va.  22217 
ATTN:  Code  715  (Tech  Lib) 


DEPARTMENT  OF  THE  AIR  FORCE 


I  Air  Force  Institute  of  Technology 

Air  University 

Wright-Patterson  AFB,  Ohio  45433 
ATTN:  Library 


1 

1 

1 


1 


1 


1 


t 

\ 


> 

i 

t 


? 

i 

j 


1 


1 


1 


t 

I 

\ 


110 


Address 


No .  of 
Copies 


WAMnm  or  the  mr  force  (o.,ti 


Air  Force  We.pon.  Ubor.tory,  AFSC 
“j^nd  AFB,  N.  Hex  87117 
ATTN:  DE  N.  Plamondon 
DES  R.  Jolley 
DES  J .  Thomas 
SUL 


Deputy  Chief  of  Staff 
Research,  Development,  &  acq. 
Department  of  the  Air  Force 
Washington,  D.  C.  20330 
ATTN:  AFRDQSM 


Norton  AFB,  Calif.  92409 
ATTN:  HNNXH  <D.  Gage) 

WRNXH  (N.  Del  Vecchio) 


Vela  Seisaological  Center 
312  Montgomery  Street 
Alexandria,  Va.  22314 
ATTN:  G.  Ullrich 
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DEPARTMENT  OF  energy 


Nuclear  Regulatory  Commission 

VwhiMtT*  iiCen,in*  Re*»i«tions 
Washington,  D.  C.  20545 

ATTN:  Site  Analysis  Br  (L.  Heller) 

Uwrence  Livermore  Laboratory 

P-  0.  Box  808 

Livermore,  Calif.  94550 

ATTN:  Technical  Information  Dept.,  Library 

Oak  Ridge  National  Laboratory 
Nuclear  Division  7 

X-10  Lab  Records  Division 
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Address 


No.  of 
Copies 


DEPARTMENT  OF  ENERGY  (Continued) 

Ssndis  Laboratories 
P.  0.  Box  5800 
Albuquerque,  N.  Hex.  87115 

ATTN:  Div  5531  (W.  Brown)  1 

Div  5532  (J.  Lipkin)  1 

Div  5621  (M.  Hightower)  1 

Div  5533  (A.  Chabai)  1 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp. 

P.  0.  Box  92957 

Los  Angeles,  Calif.  9009 

ATTN:  Technical  Information  *n  r *ices  1 

Agbabian  Associates 
250  N.  Nash  Street 
El  Segundo,  Calif.  90245 

ATTN:  M.  Agbabian  » 

Applied  Theory,  Inc. 

1010  Westwood  Blvd 

Los  Angeles,  Calif.  90024 

ATTN:  J.  Trulio  > 

Boeing  Co. 

P.  0.  Box  3707 
Seattle,  Wash.  98124 

ATTN:  Aerospace  Library  1 

California  Research  &  Technology,  Inc. 

6269  Variel  Avenue 
Woodland  Hills,  Calif.  91364 

ATTN:  Library  1 

S.  Shuster  1 

K.  Kreyenhagen  1 

California  Research  &  Technology,  Inc. 

4049  First  Street 
Livermore,  Calif.  94550 

ATTN:  D.  Orphai  I 
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Address 


No.  of 
Copies 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Civil  Systems,  Inc. 

RFD  1 

South  Royalton,  Vt.  05068 

ATTN:  S.  Blouin  1 

Civil  Systems,  Inc. 

P.  0.  Box  2083 
Midland,  Tex.  79702 

ATTN:  S.  Melzer  1 

Civil  Cvstems,  Inc. 

2201  San  Pedro,  NE 
Albuquerque,  N.  Mex.  87110 

ATTN:  J.  Bratton  1 

Eric  H.  Wang 

Civil  Engineering  Research  Facility 


University  of  New  Mexico 
University  Station 
P.  0.  Box  25 

Albuquerque,  N.  Mex.  87131 

ATTN:  J.  Shinn  1 

J.  Kovarna  1 

P.  Lodde  1 

Fugro  National,  Inc. 

P.  0.  Box  17145 

Long  Beach,  Calif.  90807 

ATTN:  D.  Anderson  1 

General  Electric  Co. -Tempo 

816  State  Street  (P.  0.  Drawer  QQ) 

Santa  Barbara,  Calif.  93102 

ATTN:  DASIAC  1 

IIT  Research  Institute 
10  W.  35th  Street 
Chicago,  Ill.  60616 

ATTN:  Documents  Library  1 
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Address 


No.  of 
Copies 


DEPARTMENT  OF  DEFENSE  C^fRACTORS  (Continued) 


University  of  Illinois 
Consulting  Services 
B106A  Civil  Engineering  Bldg. 
Urbans,  Ill.  61801 
ATTN:  N.  Newmark 

J.  H.  Wiggins  Co.,  Inc. 

1650  S.  Pacific  Coast  Highway 
Redondo  Beach,  Calif.  90277 
ATTN:  J.  Collins 

Merritt  Cases,  Inc. 

P.  0.  Box  1206 
Redlands,  Calif.  92373 
ATTN:  library 

Physics  International  Co. 

2700  Merced  Street 
San  Leandro,  Calif.  94577 
ATTN:  Technical  Library 
J.  Thomsen 

Pacifica  Technology 
P.  0.  Box  148 
Del  Mar,  Calif.  92014 
ATTN:  Technical  Library 

R&D  Associates 
P.  0.  Box  9695 

Marina  Del  Rey,  Calif.  90291 
ATTN:  R.  Port 

Technical  Information  Center 

R&D  Associates 
1401  Wilson  Blvd.,  Suite  500 
Arlington,  Va.  22209 
ATTN:  N.  Cooper 

Science  Applications,  Inc. 

P.  0.  Box  2351 
U  Jolla,  Calif.  92038 
ATTN:  Technical  Library 
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Weidlinger  Assoc.,  Consulting  Engineers 
3000  Sand  Hill  Road 
Henlo  Park,  Calif.  94025 

ATTN:  J.  Isenberg  1 


115 


In  accordance  with  Uttar  fro*  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Baladi,  I'.eorge  Youssef 

An  effective  stress  model  for  ground  motion  calculations  / 
hy  C.  V.  Baladi.  Vicksburg.  Miss.  :  II.  S.  Waterways  Experi¬ 
ment  Station  ;  Springfield.  V*.  :  available  from  National 
Technical  Information  Service.  1979. 

106  p.  :  ill.  ;  27  c*.  (Technical  report  -  U.  S.  Amy 
Engineer  Waterways  Experiment  Station  ;  Rl.-'P-'l 

Prepared  for  Defense  Nuclear  Agencv,  Washington.  D.  I'., 
under  DMA  Subtask  SB209.  Work  Unit  44. 

References:  p.  R7-89. 

1.  Constitutive  model  * .  2.  Ground  motion.  A,  Ground  shock. 

4.  Stress-strain  relations.  1.  United  States.  Defense 
Nuclear  Agency.  II .  Series:  United  States.  Waterways  Ex¬ 
periment  Station.  Vicksburg.  Miss.  Technical  report  ; 
Sl-'9-7. 

TA'.N.VI  no.  SI  ..'9 -7 


